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Abstract 

o : 

We describe the evolution of the quantity of parasites in a population of cells which 
divide in continuous-time. The quantity of parasites in a cell follows a Feller diffusion, 
which is splitted randomly between the two daughter cells when a division occurs. The 
cell division rate may depend on the quantity of parasites inside the cell and we are 
interested in the cases of constant or monotone division rate. We first determine the 
asymptotic behavior of the quantity of parasites in a cell line, which follows a Feller 
diffusion with multiplicative jumps. We then consider the evolution of the infection of 
Pj^ ■ the cell population and give criteria to determine whether the proportion of infected 

0^ cells goes to zero (recovery) or if a positive proportion of cells becomes largely infected 

(proliferation of parasites inside the cells). 
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■ 1 Introduction and main results 

C — ■ , We consider a continuous time model for dividing cells infected by parasites. We assume that 

parasites proliferate in the cells and that their lifclcngths are much less than the one of the 
cell. Informally, the quantity of parasites (X t : t > 0) in a cell evolves as a Feller diffusion 
|Fel71j . The cells divide in continuous time with a rate r(x) which may depend on the quantity 
of parasites x that they contain. When a cell divides, a random fraction G of the parasites 
goes in the first daughter cell and a fraction (1 — 9) in the second one. More generally, instead 
of parasite infection, one could be interested in some biological content which grows in the 
cells and is shared randomly when the cells divide (for example proteins, nutriments, energy 
or extrachromosomal rDNA circles in yeast |SG97j ). 

Let us give some details about the biological motivations. For the sharing of parasites, we 
are inspired by experiments conducted in Tamara's Laboratory where bacteria E-Coli have 
been infected with bacteriophage lysogens [SMPT05]. They show that a very infected cell 
often gives birth to a very infected and a lowly infected daughter cells. Thus we are interested 
in taking into account unequal sharing and we do not make restrictive assumptions about 0. 
The model we consider here is a continuous version of Kimmel's multilevel model for plasmids 
}Kim97j . In the latter model, the cells divide in continuous time at a constant rate and the 
number of parasites is a discrete quantity which is fixed at the birth of the cell: the para- 
sites reproduce 'only when the cells divides'. Moreover the sharing of parasites is symmetric. 
In |Ban08| . a discrete time model where the sharing of the parasites may be asymmetric is 
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considered. We refer to [AK 98, BP94] for general discrete time models with independent and 
identically distributed values for the offspring and to [Guy07j [DM 08 , BS GP08] for asymmetric 
models motivated by cellular aging. Here both the quantity of parasites and the time are 
continuous. 

In addition to unequal repartition of parasites in the two daughter cells, we are also interested 
in letting the division rate of the cell r{x) depend on the quantity x of parasites in the cell. 
On the one hand the cell's division rate may decrease when the quantity of parasites increases 
since the cell's stamina is affected. For instance, in |Bon06[ IHurOll IOA02j examples where 
host fecundity is reduced by the presence of pathogens are presented. On the other hand, 
increasing rates may be found in order to fight against their proliferation. Furthermore, this 
case is relevant when there is a symbiosis between the host and its parasites: parasites then 
make the rate of division of the cell increase. In this direction, we refer to |J JBR911 IJ JBR92] . 
Note that increasing functions r are also natural if we replace parasites by nutriments or more 
generally 'energy for the cell division'. Finally, we also consider the case where the division 
rate does not depend on the quantity of parasites (r constant) since some parasites do not 
affect the evolution of their host. Moreover this is mathematically natural to consider this 
simpler case first, which will be useful for non constant cases. 

The process we study is a Markov process on Galton Watson trees and [BDMT09] give asymp- 
totic results under an ergodicity assumption, which is not fulfilled here. Mathematical ap- 
proaches for the case of non constant division rates r(x) have been considered by [HH09 
for continuous time branching diffusions, when the offspring appear at the position of their 
mother. Their method, which relies on Girsanov transformations, does not hold when nonlocal 
branching is considered. In this vein, we also refer to |KE04[ [Eng07| . 

As we will observe, even when the division rate is constant, our model behaves differently 
from the discrete time model with 'constant division' of |Ban08j and the recovery criteria are 
different. 

We refer to |Daw931 IDGL02] for the nonlocal branching supcrprocess counterpart. In the 
framework of cell division, |ES07j considers a superprocess approach for cell damages, which 
corresponds to the high density limit of small and short living cells. Notice however that asym- 
metry and nonlocal branching are lost in the continuous limit. Moreover, in view of statistical 
applications |Guy07[ IDM08] , our purpose is to stick to a discrete genealogical cell tree. In the 
experiments, the number of cells can not always be considered as large. 

Let us give now a qualitative description of our model, which is rigorously defined in 
the next section. The quantity of parasites (X t : t > 0) follows a Feller diffusion (see e.g. 
[Lam67llBm76] . jRY98| (Chapter IX) and (TW89] (Chapter IV)) 

dX t =gX t dt + ^2a 2 X t dB t (1.1) 

started at xo > 0. In this work, we study the super-critical case (g > 0) with positive variance 
(a > 0). Then the parasites survive with probability 1 — exp(— gxo/a 2 ) € (0, 1). Upon survival 
their quantity goes to infinity. This model corresponds to the high density limit of a birth 
and death process for parasites with short lives and we give details of the convergence of the 
discrete model to the continuous model studied here in the next section. 
The population of cells remains discrete and each cell divides in continuous time with a rate 
r{x) which will be here monotone, measurable and bounded on compact intervals to avoid 
degenerated situations. Let us denote by Vt the set of living cells at time t, by N t = #Vt the 
cell population size at time t and by X\ the quantity of parasites in the cell i £ Vt at time t. 
Thus: 

iev t 

When a cell containing a quantity x of parasites divides, the quantity of parasites of the 
two daughter cells received by inheritance are respectively Qx and (1 — &)x, where O is a 
random variable in [0, 1] such that P(0 = 0) = P(0 = 1) = 0. Let us note that this model can 
been seen as a random fragmentation where the mass of the fragment follows a Feller diffusion. 
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A fragment with mass x splits with rate r{x) in two fragments whose mass are respectively 
equal to Ox and (1 — &)x. Let us notice similarities with BD86, BD87 for splitting intervals 
of size x at rate r(x) = x a {a > 0), which amount to considering constant quantity of parasites 
within cells. 



We aim at determining how the infection evolves in the cell population. We stress that the 
results we give in this paper do not depend on the initial quantity of parasites xq > 0. 

First, in Section [3J we determine the asymptotic behavior of the quantity of parasites in 
a cell line (Y t : t > 0). This amounts to following the infection in a cell until it divides and 
then choose at random one of the two daughter cells. This process is a Feller diffusion X with 
multiplicative jumps O occurring at rate r(.) and we prove the following extinction criterion. 

Proposition, (i) In the case where r(.) = r is constant, 

■k If g < E(log(l/0))r, then almost surely (a.s.) there exists t > such that Y t = 0. 

* Otherwise P(limt_ > _ ) _ 00 Y t = +oo) > 0. 

(ii) In the case where r is an increasing function, 

* If there exists xq such that g < E(log(l/0))r(a;o), then a.s. there exists t > such that 
Y t = 0. 

* Ifg > E(log(l/6)) sup xeK+ r{x), then P^m^ Y t = oo) > 0. 
(Hi) In the case where r is a decreasing function, 

* If g < E(log(l/©)) inf xe R + r(x), then a.s. there exists t>0 such that Y t = 0. 
-k Otherwise, P(lim t _ i . 0O Y t = oo) > 0. 

We say that the organism recovers when the proportion of infected cells goes to zero as 
time goes to infinity. But contrarily to the discrete case |Ban08[ |Guy07| , the extinction criteria 
stated above do not provide directly recovery criteria for the organism. When the division 
rate is constant, we prove in Section [5] the following criterion for a.s. recovery: 

Theorem. We assume here that r(x) — r for every x > 0. 

(i) If g < 2E(log(l/0))r, then the organism recovers a.s. in the sense that: 

hm = 1 a.s. (1.2) 

t^+oo N t 

(ii) If g > 2E(log(l/©))r then the parasites proliferate exponentially inside the cells as soon 
as the parasites do not become extincted in the sense that 

{ limsup W 1 * - ' ~ S - > 0} = {it > 0, X t > 0} a.s. (1.3) 

for every < k < g — 2E(log(l/9))r. 

Thus, when the division rate is constant, either the organism recovers a.s. or a positive 
proportion of cells becomes arbitrarily largely infected as soon as the total parasite population 
survives. Let us make two observations about this recovery criterion. First, thanks to unequal 
sharing of parasites (i.e. E(log(l/9)) 3> 1), the organism may recover a.s. although the par- 
asite growth rate is large (g 3> 1) and the cells' division rate is low (r <c 1). Second, we can 
remark the factor 2 in this criterion, which is inherited from a bias phenomenon that is typical 
to continuous time branching process and stated in (BDMT09 ; CRW9ll lHH09j . Heuristically, 
a typical cell at a large time has divided at rate 2r along its ancestral lineage. This allows 
the following surprising situation: the organism may recover a.s. although the quantity of 
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parasites in a cell line goes to infinity with positive probability. 

The situation is very different if x t— > r(x) varies. In this case, the genealogical tree of the 
cell population depends on the quantity of parasites that are evolving in the cells and we are 
no more in the frame of Markov processes indexed by Galton- Watson trees such as [BDMT09] . 
If the division rate decreases when the quantity of parasites increases, then the situation is 
quite intuitive since the cells which become too infected will keep getting more and more 
infected whereas their proportion in the cell population tends to 0. Using the infection in a 
cell line (Y t : t > 0), we give a sufficient condition for a.s. recovery (see Section 15751) . 
The case of an increasing division rate is more interesting and difficult: increase in the cell 
division rate may prevent the parasites from proliferating in the cells. There is no longer either 
recovery or proliferation of parasites and an intermediate regime appears (see the first example 
in Section [5. 2p . This leads us to introduce moderate infection, which roughly speaking means 
that the distribution of the quantity of parasites in a positive fraction of the cells alive at time 
t remains positive and bounded. The criterion for this a.s. moderated infection is conjectured 
in Section 15.2.21 and illustrated with simulations in Section [6] 

The paper is organized as follows. In Section [21 we give more formal definitions based on 
measured valued processes and state the convergence of the discrete model to the continuous 
model. In Section [21 we give the asymptotic behavior for Feller diffusions with multiplicative 
jumps in the cases of constant or monotone division rates. This gives the asymptotic behavior 
of the infection in a cell line. We derive in Section [4] the recovery criteria for constant division 
rate rates r, while non-constant cases are handled in Section [S] 

2 Preliminaries 

2.1 Measure- valued description of the infected cell population 

To describe the tree of the cell population and label its nodes, we use the Ulam-Harris-Neveu 
notation (e.g. [Daw93[ lGa!060 . The mother cell is labeled by and when the cell of label i 
divides, the two daughters are labeled by iO and i\. We introduce the set of labels: 

3 = {0}U |J{0,1}" 1 . (2.1) 

m— 1 

For every i 6 J, there exists a unique m G N* such that i G {0,1}"' and we note \i\ = m 
the generation of individual i. For all i 1 = i\ • ■ ■ i} ni and i 2 = if ■ • ■ i 2 ri2 G 3, we define their 
concatenation i x i 2 by the label i\ ■ ■ ■il ni if • • - *m 2 - We write i 1 < i 2 when there exists « 3 G 3 
such that i 2 = i x i 3 . 

Let Vt C 3 be the set of cells alive at time t and N t = #Vt be its size. For i G V t at time t, we 
denote by X\ G M+ the quantity of parasites in the cell i. The population of cells at time t may 
be represented by the random point measure on 3 x M + , Z t (du,dx) = J2iev t ^(i,xt)(du,dx). 
We define 

Z t (dx) = ^25 xi (dx), (2.2) 

iev t 

the marginal measure of Z t (du, dx) on the state space R+. 

The space of finite measures on K + , M.p{M. + ) 1 is embedded with the topology of weak con- 
vergence. For a measure fi G A4f{^+) and a positive function /, we use the notation 
(M= /) = J R+ / d^- 

The evolution of (|2.2D can be modelled by a stochastic differential equation (SDE) that 
is now introduced. This defines a solution in the space 0(R+, A^_f(K + )) of cadlag measure- 
valued processes. This space is embedded with the Skorokhod topology (see e.g. |JM86p . We 
follow in this the inspiration of |FM041 IEDMT09] . 
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Let (-B 1 , i £ 3) be a family of independent Brownian motions (BMs) and let Q(ds, du, di, dO) 
be a Poisson point measure (PPM) on R + x R + x 3 x [0,1] of intensity q(ds,du,di,d9) = 
ds du n(di) K(d6) independent from the BMs. We have denoted by n{di) the counting measure 
on 3 while ds and du are Lebesgue measures on R + . We denote by (J-t ■ t > 0) the canonical 
filtration associated with the BMs and the PPM. Then, for every (i, x) -> f(t,x) € Cl ,2 (R + x 
R + ,R) (bounded of class C 1 vat and C 2 in x with bounded derivatives), 

(Z t ,f) = f(0,x )+ f f {d s f(s,x) + gxd x f(s,x) + a 2 xdlJ{s,x)) Z s (dx)ds + M(+ 

JO JR + 

(2.3) 

/ / W 8 ,u<r{xi )(f(s,6Xi_) + f(s,(l-6)Xi_)-f(s,Xi_))Q(ds,du,di,d6), 

JO Jr + x3x[0,1] 

where xq is the quantity of parasite in the ancestor cell at t = and 

M t = I E V^xjd x f( S ,Xl)dBl (2.4) 

is a continuous square integrable martingale with quadratic variation: 

(M f ) t = [ [ 2a 2 xf' 2 (s,x) Z s (dx)ds. (2.5) 
Jo Jm + 

Proposition 2.1. For given initial condition xq, PPM Q and BM (£>*, i £ 3), there exists a 
unique strong solution to \2. S\) - ^7S\) . 

Proof. Note that as expected, the total quantity of parasites X t = L_ x Z t (dx) follows the 
Feller diffusion (|1.1[) with drift g and variance a 2 . We recall that for this Feller diffusion: 

Vf e M+, E(X t ) = x e 9t < +oo. (2.6) 

For existence and uniqueness, we use (|2.6[) and refer to similar computation in |CM07] and 
[Tra] Chapter 2 Section 2.2. 

Notice that for / £ C^' 2 (R^_, R), the process Af^ is a real square integrable martingale since 
for all t e R+: 

2a 2 \\f'\\ 2 0O E( f [ xZ s (dx)ds) < +oo. 
This ends the proof. ■ 



2.2 Prom the microscopic discrete model to the continuous model 

The continuous model defined above is the limit of a discrete microscopic model where the 
parasites follow a birth and death process. In the microscopic model that we consider, each cell 
hosts a discrete parasite population. The initial cell contains [nxo] parasites, where n £ N* 
is a parameter that will tend to infinity and where [x] denotes the integer part of x. The 
parasites are rewcighted by 1/n, so that the biomass in a cell has a constant size order. They 
reproduce asexually and die with the following respective individual rates: 

b n = no- 2 + b, d n = no- 2 +d (2.7) 

where 6, d > are such that b — d = g > 0. The cell population is described by the point 
measure 

Z?{dx) = Y J S x ?Adx) 
iev t 

where X™' 1 is the number of parasites renormalized by n in the cell i at time t. 
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Let Q 1 and Q 2 be two PPMs on R + x E x :=1 + x3xR + x R + and R+ x E 2 :=K + x3x K + 
with intensity measures ds n(di) du and ds n(di) du K(d8). We associate Q 1 to the cell divisions 
while Q 2 corresponds to the births and deaths of the parasites. 

Z t = $(®,[nx ]/n) 



J E 



JE; 



Q 1 (ds, di, du, dO)\ eV ^ ^ u <r(X^) (&(il,[BnX?^]/n) + \i2,X^-[$nX^]/n) \i,X^)J 

Q 2 (ds,di,du)l ieVs _ [[S (i XT A +1/n) - S^x^)) K 



i<b n X r s 



^(i,X?*-l/n) \i,X™^) J \ n X"^ <u<{b n +d n )X r ' 



(2. 



Note that other discrete models would lead to this continuous model. For example, parasites 
may have several offspring's and their sharing in the two daughter cells could follow a random 
binomial distribution. This means that we draw 6 in the distribution K{d6) and send each 
parasite in the first daughter cell with the probability 9, and else send it in the second daughter 
(see also |Kim971 lBan08| ). 

The parasite population is a continuous time birth and death process of constant rates b n and 
d n . Hence for every n £ N* , there is existence and strong uniqueness of the solution of (|2.8|) 
for a given initial condition Zg and PPMs Q 1 and Q 2 . 

Proposition 2.2. Assume that there exists an integer p > and a positive f > such that for 
all x € R+, < r(x) < f(l + x p ). Then, the sequence (Z n : n £ N*) defined in \2. 8\) converges 
in distribution in B(R + , A4f(^+)) as n — >■ +oo to the process Z defined in i2.3\) - ^75\) . 

The proof is deferred in Appendix (Section [7]). Additional regularities on the division rate 
r would be required to control the difference between the microscopic process (|2.8p and its 
approximation (|2.3|) - (|2.5p . 



3 Feller diffusion with multiplicative jumps 

We are interested in the evolution of the quantity of parasites in a cell line. This means that at 
each division, we only keep one cell and consider the quantity of parasites inside. This process 
(Yt : t > 0) follows a Feller diffusion with multiplicative jumps whose sizes are independent 
and distributed as the r.v. and which occur at rate r(x) when the process is equal to x. 
Without loss of generality, we decide to follow the quantity of parasites in the 'first branch' of 
the tree, which is constituted by individuals labeled only with 0s. More precisely, letting 

n = {0} u |J {or, 

m£N* 

this process (Y t : t > 0) is defined by: 

Y t = Xl o[t \ where i (t) e 1 D V t . 

Using (|2.3[) . it satisfies 

Y t =x + I gY s ds + f v /2a 2 Y s df3 s 
Jo Jo 

where 

d/3 t = $z(i)dBl and p(ds, du, d6) = Ibnv s _ (i)Q(ds, du, di, d6) 

iev t 

are respectively a BM and a PPM on E := M + x M + x [0, 1] with intensity ds du K(d6). In the 
sequel, we denote by (J 7 ^ : t > 0) and (J 7 ? : t > 0) the canonical filtrations associated with 



-x[0,l] 



\< r (Y s _) (1 - 0)Y s _p(ds, du, dS) (3.1) 
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the BM /3 and the PPM p respectively. 

Apart from the biological motivations considered here, one can notice that such Markov 
processes with multiplicative jumps have various applications (see for instance the TCP in 
|DGR021 IGRB041 UatsSj 'l. 

3.1 Extinction criterion when r is constant 

We first study the asymptotic behavior of the process Y when the jump rate does not depend 
on the value of the process. 

Proposition 3.1. We assume that r(.) = r is constant. 

(i) Ifg< E(log(l/9))r, then P(3i > : Y t = 0) = 1. 
Moreover if g < E(log(l/0))?\ 

3a > 0, Vx > 0, 3c> 0, P Xo (Y t > 0) < ce~ at (t > 0). (3.2) 

(ii) Ifg > E(log(l/6))r, then P(Vi > : Y t > 0) > 0. 
Furthermore, for every < K < g — E(log(l/©))r > 

{ lim e - Kt Y t = oo} = {Vi : Y t > 0} a.s. (3.3) 

To guess this extinction criterion one can observe that without division the parasite pop- 
ulation follows a Feller diffusion with ¥,(X t ) = xocxp(gt). With the multiplicative jumps 
corresponding to the cell divisions, we obtain: 

M 

Y t w x e gt Y[ ~ x cxp (gt +jV ( x E(log(0))) (t -> oo) 
i=i 

where Aft is a Poisson r.v. of parameter rt and where the 0/s are i.i.d. r.v. with distribution 
K{d9) and independent of Aft- 

More rigorously, to prove the proposition above, we compute the Laplace transform of the 
jump-diffusion process (|3.1[) . Let us introduce the following rescaled process corrected with 
its drift and jumps: 

Y t =Y t e K \ 

with 

M 

3= ' 

Lemma 3.2. The process (Y t : t > 0) is a continuous local martingale and for all t, A, Xq > 0, 



K t = -gt - l u < r \og(9)p(ds, du, d9) = -gt - V log^ 

Jo Jr+x[o,i] - =1 



Proof. Using that for every i G R+, < Y t < X t and (|2.6p . all the stochastic integrals that we 
are writing are well defined as local martingales. Using Ito's formula with jumps {e.g. Ikeda 
Watanabe Th.5.1 on p.67 [IW89] ) : 

Y t =x + [ e Ks \gY s ds + ^2a 2 Y s d(3 s ] - [ gY s e K °ds 
Jo L o 

+ / / (Y s e K ° - Y s _ e K '- ) K u < r p(ds, du, d6) 

Jo il + x[0,l] V ' 

=x + f e K ^2^Y s dp s + f f F s „ (fcT log(e) - l) % t < r p(ds, du, dB) 

Jo Jo 7r + x[o,i] v ' 



--xo 



j e K W2o- 2 Y s d[3 s 
Jo 
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Then (Y t : t > 0) is a continuous local martingale which satisfies: 

ft 



Yt=x a + e K ^W^ 2 Y s d(3 s . (3.6) 



o 



Let us now work conditionally on J-^ . 

Using Ito's formula for a function u(t, y) which is diffcrcntiablc by parts with respect to t and 
infinitely diffcrcntiable with respect to y, we get 



u(t,Y t ) = u(0,x o )+ 



^(s,Y s )e K °/ 2 V2^Ysdf3 s 
dy 



Following Ikeda and Watanabe }IW89] (8.10) on p. 236, we choose u(t, y) so that the finite 
variation part equals zero. More precisely, letting to > and 

u(t,2/):=exp( t \ ■ — | (0 < * < t ), (3.7) 



a 2 X J; e K -ds + l j 



so that 



«(*,«) = «((),*„)- / ^ + i exp ^- - . ^ + ] df3 s . 



The process (tt(t, Y t ) : < t < to) is a local martingale bounded by 1 and thus a real 
martingale. We deduce that 

E X0 (u(t , Y t0 ) | J£) - E K0 («(0, F ) | J£) , 

which gives 

E X0 (exp(-AF t0 ) | 7*) = exp (- Ax ° ) . (3.8) 

\ a A J e^ s ds + 1/ 

This provides (|3.5|) by taking the expectation. ■ 

We can now prove Proposition 13 . 1 1 in the case of constant division rate. 
Proof of Provosition HOI The map 1 1— ^ J * exp(J4T s )ds is an a.s. increasing function and thus 

/•+00 

Zoo := / exp(K s )ds G M+ U {+00} 
Jo 

is well defined. Using dominated convergence, the r.h.s. of (|3.5[) converges to 

E, (exp(-A.T /( ( 7 2 AZ oo + l)). 

Using Lemma [XU the process (Y t : t > 0) converges in distribution as t — > +00 to Y"oo whose 
distribution is specified by 

^o(e-^)=E, (ex P (-^ TT )). (3.9) 

Recalling from (|3.6[) that (Yj : t > 0) is a positive local martingale, wc obtain by Jensen's 
inequality that (exp(— Y t ) : t > 0) is a positive sub-martingalc bounded by 1. From this, we 
deduce that the convergence towards Yao, which is possibly infinite, also holds a.s. 
Letting A — ¥ +00, we get by bounded convergence: 

p -ff» = 0) =Aft»^ ( cxp (-^Tl)) = E *° ( exp (-^)) ' (3 ' 10) 
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From |Ber96j Corollary 2 p. 190, we know that the Levy process (K t : t > 0) defined in 
tends to +oo (resp. — oo, resp. it oscillates) when E(Ki) = E(log(l/G))r — g is positive (resp. 
negative, resp. zero). Thus, we split the proof into three cases. 

If g > E(log(l/6))r, we choose k > such that g-K > E(log(l/6))r. Then (K t +nt : t > 0) 
is a Levy process such that E(Ai + k) < 0, so 

lim K t + at = —oo a.s. 

t— >-+oo 

We define B := sup t>0 {A t + Kt} < oo a.s. so that K t < —Kt + B a.s. This ensures that 
Zoo < oo a.s. and using p.lOp . we get: 

P(Foo=0)<l. (3.11) 
As Y t = Y t exp(— K t ) and lim t _j. +00 exp(— K t ) = +oo, we have: 



P Nt > : Y t > 0) > P lim Y t = +oo > P lim Y t > > 0. 

Furthermore, for every < K < g — E(log(l/9))r, —K t — Kt — ^ +oo as t — > oo. So A = 
{lim t _5. +00 K t + Kt — +00} G J-?, happens a.s. and 

lim e- Kt Y t = +00 I ) > P f lim F t > | T 1 ^ ) > a.s. 

by using (|3.8[) as in (|3.1ip . We can now prove that {lim t _j. +00 e~ Kt Y t — +00} — {lim t ^ +00 Y t = 
+00} a.s. Let TV > 0. Conditionally on the event {Y t — > 00}, the stopping time T/y = inf{i > 
: Y t > TV} is finite for every TV. Since Pi (lim t ^ +00 e~ Kt Y t = +00 | J 7 ^) > a.s., and since 
the process (Y t : t > 0) satisfies the branching property conditionally to J-^ , this ensures that 



N 



lim e~ Kt Y t < +00 \TfL \ < Pi lim e - Kt F t < +00 | ) ^0 

t— > + C30 / \ t— > + 00 



as TV — >• +00. Thus conditionally on J 7 ^ and {Y t — > +00}, e~ Kt Y t — > +00 a.s. Finally the 
fact that {lim t ^ +oc Y t — +00} — {V< : Y t > 0} a.s. is a classical consequence of the Markov 
property using that is an absorbing state. 

If g = E(log(l/8))r, the Levy process (K t : t > 0) oscillates a.s.: 

limsup Kt = +00, liminf K t = — 00. (3-12) 

Then, for every k > 0, the stopping time :— inf{t > : K t > k} is finite a.s. This ensures 
that almost surely 



Zoo > / 1 e Kt dt>e k [ ' e Kt ~ KT * dt > e k exp ( inf {TV t - K T J). (3.13) 
Jt,. Jt.. te[T fc ,T fc +i] 



r T k + l r T k + l 

As (inf 4g [ Tfcj j' fc+ i]{iV t — Kx k } ■ k > 1) arc identically distributed finite r.v., we have 



lim sup e fc exp ( inf {iV t — AV i) = +00 a.s. 

kew y te[T k ,T k +i] 1 k " 

Since the l.h.s. of (|3.13[) does not depend on fc, letting k — > 00 ensures that Z^, = +00 a.s. 
and l|3.10|) gives: 

Yoo = a.s. (3.14) 

Our purpose is now to prove that (Y t : t > 0) reaches in finite time a.s. Let us define the 
following stopping times, which are finite from (|3.12p and (|3.14[) : 

Ti := inf {i > : Y t < 1}, r l+1 := mi{t > n + 1 : Y t < 1} (i > 1). (3.15) 
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Introducing a := ini^o ^jP^Yi = 0)} > we have: 

j 

P(Y Tj > 0) = P(Y T1 > 0) JJP(F Tj > | Y Ti _, > 0) 

i=2 
3 

< P(Y T1 > 0) JJP(y Ti _ 1+1 > | Y Ti _, > 0) < P(Y T1 > 0)(1 - a y-\(3.i6) 

i=2 

which tends to as j — > oo. Then, P(3n £ N : Y Tri = 0) = 1, which is the desired result. 

Finally if g < E(log(l/9))r, we choose k > so that g + n < E(log(l/9))r. Then K t -nt 
+oo a.s. Proceeding as in the case g > E(log(l/0))r with B := mi t ^R + {K t — nt] > — oo a.s., 
we obtain K t > nt + B a.s. This implies that = +oc a.s. Using (|3.10[) . we get Y x = a.s. 
Moreover exp(— K t ) — > as t — > oo a.s., so 

lim Y t = lim Y t e- Kt = a.s. (3.17) 

t— >+oc t— >+oo 

Using again the stopping times (|3.15[) . we get that Y t reaches in finite time a.s. 
Let us now prove (|3.2[) . Formula (|3.5[) yields : 

P, o (F t = 0)= lim Efexpf- - . A *° ; \ \ = E f cxp I H 



A^+oo ^ ^ a 2 \f Q e K sds + l)) \ \ a 2 f Q e K °d S/ 
As the process {K t '■ t > 0) has no negative jumps and drift —g, we have 



for every i > 1 and 



inf if„ > — g a.s. 

«e[i— 



e*"ds > / e^ds > e^- 1 " 9 a.s. 
./t-i 



Moreover for all x > and a G [0,1], 1 — e 1 < min(l,a;) < min(l,a; Q ) < x a . This gives for 
every t > 0, 

P Ko (y t >0) = E(l-cxp(-x a- 2 e- Kt - 1+9 )) 
< (x 0f r- 2 e s ) a E (e"^- 1 ) 

using the Levy-Khintchine formula where <j> is the Laplace exponent of (K t ) t >o (see |Ber96j ): 
<f>{a) := -5a + rE(l - e - Ql °s( 1 /e)) ( a > )_ 

Adding that 0(0) = and <f/(0) = rE(log(l/9)) - g > 0, there exists a G (0,1] such that 
4>(a) > and: 

P ao (F t > 0) - P X0 (Y t > 0) < [(a:o<r- 2 e ff ) a e*W]e- t *W > 
which completes the proof. ■ 



3.2 Extinction criteria with monotone division rate 

We give here the extinction criteria of the process (Y t : t > 0) describing the quantity 
of parasites in a cell line when the jump rate r is monotone. For the proof, we use coupling 
arguments to compare this process with the case of constant division rate. 
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We begin with the case where r is an increasing function which means that the more 
parasites the cells contain, the faster they divide. This case is relevant when the cell division 
rate is increased to get ride of the parasites or when there is a symbiosis between parasites 
and cells. The asymptotic behavior of Y depends on the maximum division rate 

r* := sup r(x). (3.18) 

Proposition 3.3. We assume that r is an increasing function. 

(i) If there exists x\ > such that g < E(log(l/0))r(xi), then 

f(3 t > 0, Y t = o) = 1. 

(ii) If g > E(log(l/0))r*, then P(Vi > : Y t > 0) > 0. 
Furthermore, for every < K < g — E(log(l/0))r*, we have a.s. 

{ hm e~ Kt Y t = oo} = {it > : Y t > 0}. 

t— 5- + 0O 

Let us note that the case g > E(log(l/0))r(x) for every x > and g = E(log(l/0))r* 
remains open. The expected result is a.s. extinction but this may depend on the speed of 
convergence of r{x) to r* as x — > oo. 

Proof. Heuristically, if g < E(log(l/0))r(xi), as soon as Y > x±, the division rate is larger 
than r{xi) and Proposition 13. II ensures that the process is pushed back to x\. Eventually, it 
reaches zero. 

We give the proof of (i) using a coupling argument. Let us define Y as: 

~ /•* ~ r t+1 i ~ — r t+1 r ~ 

Y t =Y 1 + gY s ds+ V 2^Y s ^dB a s / l u < r(xi) 6Y s __ lP {ds, du, dd) 

Jo Ji Ji Ju + x[o,i] 

with initial condition Yi, with the same BM f3 and PPM p as Y shifted from 1, and with the 
constant division rate r(x) = r(x\). 

Thus, Y is a Feller diffusion with multiplicative jumps given by and constant division rate 
r(xi). Proposition 13.11 ensures that it becomes cxtinctcd in finite time. 
Moreover, denoting by 

ri := inf{< > 1 : Y t < xi} G [0, oo], 
the definition of Y ensures that 

Y 1+t <Y t a.s., (0<*<n), 

since Y undergoes less jumps than Yi+.. As Y becomes cxtinctcd in finite time, T\ < oo a.s. 
Similarly, the following stopping times: 

t :=0, Tj+i := inf {t > n + 1 : Y t < xt} (3.19) 

are finite a.s. for i e N. Proceeding as in (|3.16[l . we obtain that ¥(3i € N : Y Ti = 0) = 1. This 
ends the proof of (i). 

We now prove (ii) and we assume that g > E(log(l/0))r*. We define (Y t : t > 0) with 
the same BM and PPM as (Y t : t > 0) except that the indicator in (|3.1j) is replace by ll u < r * . 
The definition of r* and this pathwisc construction ensure that for every t > 0, Y t > Yt 
a.s. Moreover Y is a Feller diffusion with multiplicative jumps with constant rate r* and 
Proposition 13.11 (ii) states that Y grows geometrically with positive probability. Thus the 
same holds for Y. Combining the Markov property with Proposition 13.11 (ii) ensures that 
{limt-^-i-oo exp(— Kt)Y t } = {\/t > : Y t > 0} a.s. since we have the analogous result for the 
coupling process Y. ■ 
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We now consider the case where the more parasites there are, the less the cells divide. This 
is natural if the parasites make the cells ill or use their nutriments. The asymptotic behavior 
now depends on 

r* = inf r(x). (3.20) 

x>0 

Proposition 3.4. We assume that r is a decreasing function. 

(i) Ifg < E(log(l/0))r„ then v(3t > 0, Y t = o) = 1. 

(ii) Else, P(Vi > : Y t > 0) > and for every < n < E(log(l/6))r* - g, we have a.s. 

{ lim e~ Kt Y t = oo} = {Vt > : Y t > 0}). 

t— >oo 

Proof. Let us define the process (Y t : t > 0) with the same BM and PPM as (Y t : t > 0) but 
we replace the indicator in (|3.1[) by ]l u < r , . Then, for every t > 0, Y t > Y t a.s. and Y is a 
Feller diffusion with multiplicative jumps and constant jump rate r*. If g < E(log(l/8))r*, 
Proposition 13.11 ensures that Y becomes extincted a.s., which entails (i). 

For (ii), and let us consider x\ such that g > E(log(l/0))r(xi). We define the process 
(Y t :t > 0) with the same BM and PPM as (Y t : t > 0) but we replace the indicator in ([33]) 
by i u < r ( Kl ), so that it divides with the constant rate r(x\). Using the Markov property and 
Proposition ^. 11 we get 

P( lim Y t = oo) >P(Yi >xi + l)¥ Xl+l { lim Y t = oo) 

t— too t^-oo 

>¥(Y 1 > Xl + l)P* 1+ i( lim Y t = oo; Vt > 0, Y t > Xl ) > 0, 

t— >oo 

which gives the first part of (ii). The second part comes from the Markov property and this 
coupling argument as for the two previous propositions. ■ 



4 Recovery criterion for constant division rate 

In this section, we want to determine how the infection evolves in the cell population. More 
precisely, we are interested in the asymptotic proportions of cells which contain a given quantity 
of parasites. 

The questions we focus on do not need spatial structure on the cell population so without loss 
of generality we assume by now that O is symmetric in distribution with respect 1/2: 

9 = 1- 6. 

In the case of a constant division rate r, (N t : t > 0) is a Yule process and E(7V t ) = e rt . For 
the recovery criterion of the organism, we are interested in the asymptotic behavior of 

Z t {dx)/N t . 

But, as usual for branching Markov process [BDMT09, ESQ7], it is more convenient to consider 
the following renormalization : 

lt {dx) := Z t {dx)/E(N t ) = e - rt Z t {dx). 

Actually, we just need here the expectation of this quantity, whose evolution is given by the 
following result 

Lemma 4.1. The family of probability measures (E(7 t ) : t > 0) is the unique solution of the 
following equation in [y t : t > 0) for (/ : (t, x) M- ft{x)) S C^' 2 (M+ x R + ,R) and t £ R + : 

(y u ft) = fo(x ) + 2r [ [ f [f s (6x) - f s (x)] K{d0) v s (dx) ds 

JO JR + JO 

+ [[ (d s f s (x)+gxd x f s (x) + o- 2 xd 2 xx f(x))v s (dx)ds. (4.1) 

JO JR 
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This result ca be derived directly from the measure-valued equation (|2.3p and we give the 
proof below. We can then interpret E(7t) as the marginal of an auxiliary process (£t : t > 0): 

V/eC 5 2 (M + ,K), VteK + , (E( 7t ) ! /) = e- rt E(^/(Xl))=E(/(6)), (4.2) 

iev t 

where : t > 0) is a Feller diffusion which jumps with rate 2r from a; to Qx. More precisely, 
it is defined for t > by 

tt = *o+ [ g€,d8+ [ y/2^T.dW a + [ f iu<2r[f(0£ s J-f(£s-)}N(ds,du,d6) 
Jo Jo Jo Jmx[0,1] 

(4.3) 

where N(ds, du, dO) is a PPM of intensity du dO and W, a standard BM independent from N . 

This result is generalized for Markov processes indexed by Galton- Watson trees in [BDMT09] , 
with a different approach which leads to pathwise representation. The auxiliary process jumps 
with rate 2r whereas the cell divides with rate r. This bias phenomenon is classical and have 
been obtained in [CRW91, IHH09j with different approaches. It corresponds to the fact that 
the faster the cells divide, the more descendants they have at time t. That is why the ancestral 
lineages from typical individual at time t have an accelerated rate of division 2r. 

Using dominated convergence in the l.h.s. and r.h.s. of (|4.2[) . we show that (|4.2[) also holds 
for f{x) = IfcX). By Proposition 13.11 we can now determine the evolution of parasites in the 
cell population and prove the recovery criterion, when the division rate is constant. 

Theorem 4.2. (i) If g < 2rE(log(l/0)), then the organism recovers a.s. in the sense that: 

y #{* g Vt ■ Xj = 0} , 

lim = 1, a.s. 

t->+oo N t 

(ii) If g > 2rE(log(l/0)) then the parasites proliferate inside the cells as soon as the parasites 
do not become extinct in the sense that 

{ l im sup #{» £ V t .XI > e^} > q} = {\ft > : X t > 0} a.s. (4.4) 

t-too Nt 

for every < k < g — 2rE(log(l/0)). The probability of these events equals l — exp(—gxo/a 2 ). 

The factor 2 in the criterion comes from the auxiliary process and 'increases recovery' in the 
sense that the quantity of parasites in a cell line Y may go to infinity with positive probability 
whereas the organism recovers a.s. Note that there is a zero-one law: either the cells recover 
or parasites proliferate inside a positive proportion of cells. This dichotomy may fail when r 
will be an increasing function of the quantity of parasites (Section 13. 2[) . 

Proof of Lemma\Jl\ Let t e R+ and (/ (s, x) i-> f s (x)) e Cl' 2 (R + X K+,K). Using (£3} with 
(s, x) i-> f s (x)e~ rs entails: 

{it, ft) = fo(x )+ [ f (gxd x f s (x) + a 2 xd 2 xx f s (x) ~ rf s (x) + d s f s (x)) e - rs Z s (dx) ds + M t f 
Jo Jr 

+ 11 %ev s lu<r [fs(0Xl_) + / s ((l - 6)Xl_) - f s {Xl_)] e- rs Q(ds, du, di, dO) 

Jo J3xRx[0,l] 

where M{ is a continuous square integrable martingale started at 0. Taking the expectation: 

<E( 7t ),/ t ) =f {x ) + f f (gxd x f s (x)+a 2 xd 2 xx f s (x)+d s f s (x)) E( ls )(dx)ds 
Jo Jm. 

+ [[ f 2r{f s (6x)- f s (x)]K(d6)E( ls )(dx)ds, (4.5) 
Jo Jwt+ Jo 
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by using the symmetry of K(d&). 

Let us prove that there is a unique solution to (|4.ip . Let (y\ : t > 0) and (y\ : t > 0) be 
two solutions. Recall (e.g. |Rac91j 'l that the total variation distance between v\ and v\ is 

\\v\-v\\\ TV = sup (4.6) 

06C(R + ,R) 
IMI=o<l 

Let t G R + and (p G Cf(K+,R) with H^H^ < 1. We denote by (P a : s > 0) the semi-group 
associated with the Feller diffusion (jl.ip started at a: 6 R+: P s ip(x) = K x (<p(X s )). Notice that 
ll-Ff-s^lloo < IMloo < 1- Then using f|4.5[) with f s (x) = P t - S f{ x )i the second term equals 0: 

\{i$-$,<p)\ = 2rJ^ jf ^ (p t _ sV3 (te) - ^^(^^(^(i/i - y , 2 )(&)& 

<4r f \\v] - v 2 s \\ TV ds. (4.7) 
Jo 



Since C^(R+,R) is dense in Cf,(R+,R) for the bounded pointwise topology, taking the supre- 
mum in the l.h.s. implies that: \\v} — v^\\tv j^\\ v l~ v 1\\Tvds. Gronwall's lemma implies 
that \\u] - v%\\tv = 0. ■ 

For the proof of the theorem, the following Lemma will be used to obtain the almost sure 
convergence from the convergence in probability: 

Lemma 4.3. Let V be a denumerable subset and (N t (i) : t > 0) be i.i.d. processes distributed 
as the Yule process (N t : t > 0) for i G V . Then there exists a nonnegative nonincreasing 
function G on R+ such that G(y) — > as y — > oo and for all I, J finite subsets ofV and x > 0: 

^ m >*)<Gm*). (4-8) 



»(sup^ 



ieI N t (i) ~ J ~ 
Proof. We introduce for every i G V, 

M(i) := svcpN t (i)e~ rt and m(i) = inf N t (i)e 



t>o 



t>o 



As (N t : t > 0) is a Yule process, (M(i) : i 6 V) and (m(i) : i G V) are both finite positive 
i.i.d. r.v.'s with finite expectation. Moreover 



Ei e ^t(i) " E i6 /m(<) " #/ #J E< 6/ m(i)' 
and the result follows by defining for y > 

G(y) = sup (p( E ^f W # 7 > ) . j j c y. #7j # j < oo]. 
Indeed, by the law of large numbers, the sequence 



is uniformly tight. So G(y) — > as y — > oo. ■ 

Proof of Theorem \4-2\ Let us start with the proof of (i). We can apply (|4.2p with f(x) = i x >o 
since this function is the increasing limit of function which belong to / G C^(R+,R). Then 
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where (£ t : t > 0) is defined in (|4.3p . By Proposition 13.11 under the assumption of (i), 
(£t : t > 0) becomes extincted in finite time a.s. so that, (jt, ^x>o) converges in L 1 and hence 
in probability to 0. As K(N t )/N t tends in probability to 1/W where W is an exponential r.v. 
of parameter 1 (e.g. Athrcya and Ncy AN70 Chap. Ill Sect. 4), then 

^t:I t '>0} ^ in probability. (4.10) 

To get that the convergence holds a.s. and complete the proof of (i), let us denote by V t * = 
{i G Vt : X'l > 0} the set of infected cells and by N* = #V* its cardinal. Denoting by Vt iS (i) 
the set of cell alive at time t + s whose ancestor at time t is the cell i £ Vt, we have 

where (Vt, s {i) : s > 0) are i.i.d. for i G Vt. As Nt — > oo and N£ /N t — > as t — > oo, we get by 
Lemma 14.31 that 

S«Ey t *#*U0 n . ...... 

lim sup -= — — = I) m probability. 

t-H-oo s > ^ ieVi #Vt, s {i) 

Moreover for all e,7y > 0, there exists t > such that F(N^/N t > rj) < e. Then, choosing r/ 
small enough, we get that 

P(supA t * +s /A t+s > 2e) < v(sup ^ ieV "*^' S ; l) > 2e) < 2e. 

This gives the a.s. convergence and ends up the proof of (i). 

Let us now prove (ii). If there exists k G [0,g — 2rE(log(l/0))) such that: 

P(hmsu P * {l £ V * : *j ^ GXP( ^ )} > 0) = 0. (4.11) 

then lim t _ ) . +00 #{i G Vt : X£ > exp(Ki)}/7V t = in probability. Since N t /E(N t ) converges in 
probability to an exponential r.v. W of parameter 1, then 

hm — — - 1 ~ v " = in probability. (4.12) 

t->+oo E(iVt) 



Moreover 



#{i G V t : X* > exp( K i)} < N t 



E(N t ) -E(JVt)' 



which is bounded in L 2 . Then (74, fl 2; > CX p( K f)) = #{« G Vt ■ X\ > exp(Ki)}/E(A t ) is uniformly 
integrable and the convergence in probability of (|4.12j) implies the L 1 convergence. Thus, 

lim P(£ t > exp(nt)) = lim E (( 7t , H x > exp («i))) = 0, 

which is in contradiction with Proposition 13.11 (ii). Then (|4.11[) does not hold and for every 
KG[0, 5 -2rE(log(l/e))), 

Wl #{t G V t : Af > exp(Kt)} N 

P( lim sup — -f-= — > 0) > 0. 

Note that P(Vt > : X t > 0) > and by a zero one law argument, we prove now that 

{lim sup * {l £ Vt ■ X J T ~ CXp(Kt)l > 0} = {Vt > 0,A t > 0} a.s. (4.13) 

t-voo Nt 
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In that view, let us define 

V? = {teV t : x\ > 1} 

the set of cells at time t whose quantity of parasites is more than 1. Let us note that g > 
2E(log(l/F))r > r so the exponential growth of parasites is larger than the exponential growth 
of the number of cells. Then conditionally on the survival of parasites {it > Q,X t > 0}, the 
number of cells whose quantity of parasites is more than 1 can not remain bounded: 

P(limsup#^ 1 = oo) = 1. (4.14) 

t— ► oo 

Indeed if limsup t _ ) . 00 #V t 1 < oo and X t grows exponentially with rate g, then for every I > 0, 
we can find some cell Cj which contains more than / parasites. Moreover, denoting by iVj(t) 
the number of cells at time t which are issued from the same given original cell with / parasites 
and whose quantity of parasites at time t is more than 1, we have 

lim sup Nj(t) = oo. 

This gives (|4.14[) . Then the stopping time 

T n := inf{< > : #Vt > n } 

is finite a.s. Denoting by Vt(j) the set of cells alive at time t when the root of the tree is taken 
in j, we have 

#{» € V Tn+ t ■ X* Tn+t > exp(«t)} > ^ #V t (j) #{j e V t (j) : Xj(j) > expjnt)} 
N Tn+ t ~ .frf N Tn+t #Vt(j) 

Then letting t — > oo in this inequality and noting that #Vr„ > n gives 

msup ^e^x^exp^ = Ql ^ Q:Xt>0) 

t-yoo N t 

< P l( limsu P #{* g Vt ; * > exp(.t)} = B 

Letting n — > oo and recalling that the first part of the proof ensures that P X0 (limsup t _ ) . oo #{i £ 
V t : X\ > exp(nt)}/N t = 0) < 1 leads to 

PQimsup *^ £ V > : X l > g^Ml = Q | yt > : X t > 0) = 0. 

This ensures that (|4.13j) holds. ■ 



5 Evolution of the infection with variable division rate 

We now turn to the case of a variable division rate r(x), meaning that the cell division depends 
on its infection. We are still interested in the proportions of cells with a given number of 
parasites. The consideration of Z t (dx)/E(N t ) as in the previous section is not useful any 
longer and does not give simplification. Thus, we first give an SDE for the evolution of 

IH = Zt(dx)/N t , (5.1) 

Using this equation and the asymptotic behavior of the quantity of parasites in a cell line, we 
give then some asymptotic results and a conjecture for the case where r is monotone. 
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5.1 Evolution of the proportions Z t (dx)/N t 

We begin with writing down an SDE for the evolution equation for fi t when t > 0. 
Proposition 5.1. For every f e C 2 (M + ,R), 



( Mt , /> =/(*o) + / / (gxf'(x) + xa 2 f"(x)) n a {dx) ds + + M t 2 '^ 



2r(x) 



N s + l 



f(8x)K(dO) - f{x) 
f(y)K s {Hs,dy) - f(x) 



lo N + 1 

where K s (fi s , dy) is the probability measure characterized by: 
V/eC 6 (R+,R+), / f{y)K 8 {ix s ,dy) 



jj, s (dx) ds 
fi s (dx) ds, 



(5.2) 



r(y) 



If 2 f 

and where M t and M t are two martingales with quadratic variation: 

rt x 



(M 2 'f) 



o N 



iii s {dx), 2a 2 xf' 2 (x))ds 



o (^v s + i) 2 y R+ 



r(x) / {f (6x) + /((l - 9)x) - f{x) - (fi s , f)) 2 K{dd) ^(dx) ds. 

(5.3) 



Proof of Proposition \5.1l The number N t of cells alive at time t, solves the following SDE: 



N t =1 + 



^l + x3x[0,l] 



ligy 3 _ H tl <r(x* )Q(ds,du,di,d6). 



(5.4) 



From (|2.3[) and (|5.4[) . we obtain by Ito's formula for processes with jumps (e.g. Ikeda Watanabe 
Th.5.1 p.67 HUSH) that for every / G C 2 (]R + ,]R), 



(ZtJ) f( , 



where 



M, 



i,/ _ 



N 



is a continuous square integrable martingale with quadratic variation: 



{M 1J } 
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2a 2 xf'(x) 2 ^^ 



ds = 



2a 2 xf'(x)' 
N» 



[i s (dx) ds, 



and J( is the jump part: 



JR + xOx[0,l] 



ieV 3 _, u<r(X*_) 



f(Z s -,f) + f(6x) + f((l-9)x)-f(x) 



N s - + 1 

(Zs-J) 



N 



Q(ds, du, di, d9). 



17 



f f 2 f f 

Let us denote by J t = V t + M t the semi-martingale decomposition of J t . We aim at 
rewriting (|5.2[) to let the infinitesimal generator of a Markov process appear. The finite 
variation part V/ of j{ rewrites as: 



V 



r(x) 



(Z S J) jlf{9x)K{dB) f(x) (Z a ,f) 
N s + 1 N s + 1 N s + 1 iV s 



Z s (dx)ds (5-5) 



by using the symmetry of K(d9). As 

(Z a ,/> (Z s ,/)_ (Z.,/> 



iV s + l 7V S A S (A S + 1) N s + T 



we obtain: 



r(x) 



N. 



2r(x) 



N s + 1 



2 / f(6x)K(d9)-f(x)-((i s ,f} 



jjL s {dx) ds 



f(6x)K{d9) - f(x) 



fj, s (dx) ds + A, 



where the last term of (|5.6[) rewrites 



A = 



N 



( x ) ^ 1 (f(x) - (ii a ,f))fj,s(dx) 



ds 



N„ + 1 



Jy s 



N s + 1 



r ( x ){f{x) - f{y))Vs{dy)fM s (dx) ds 



(»s,r) (f -^^fix^dx) - f(y) j Ms (dy), 



by using the Fubini theorem and the fact that fi s is a probability measure. 
The bracket of the martingale part is: 



(M 2 >f) t 



4 f . iV 



1 f (gg, /) + + /((I - - fix) {Z s , /) 

A^ s + 1 N s 



(5.6) 



Kid6)Z s idx)dt 



{Ns + l) 2 7 



(/(fe) + /((l - fl)a:) - /(x) - (/i s , /)) 2 K(d6)fis{dx)ds. 



This achieves the proof. ■ 

A probabilistic interpretation of the Markov generator in (|5.2[) is in progress |BGT10j and 
described further after (I5.7I) . 



5.2 Moderate infection for increasing division rates 

We assume here that r(x) is an increasing function of the quantity of parasites. This means 
that the more the cell is infected, the faster it divides. Low infected cells divide slower and 
may even stop dividing if r(0) = 0. That's why a new regime appears here, between recovery 
and proliferation of the parasites, where a positive fraction of cells is infected but the quantity 
of parasites inside remains bounded. We then say that the infection is moderated. First, we 
provide two examples where the infection is indeed moderated: the organism does not recover 
but the parasites do not proliferate in the cells. Then, we conjecture a criterion so that the 
proportion of cells infected by more than A parasites tends to zero as A goes to infinity. We 
illustrate this conjecture with simulations in Section[6] Actually 'moderate infection' could be 
defined in several ways and we consider then the average quantity of parasites per cell. 
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5.2.1 Two simple examples 

Wc introduce two simple examples, which exhibit new behavior. First, we consider the case 
where 

6 = 1/2 a.s., r(x) = if x < 2 and r(x) = oo if x > 2. 

As soon as the quantity of parasites in a cell reaches 2, the cell divides and the quantity of 
parasites in each daughter cell is equal to one. The parasites do not proliferate in the cells 
since the quantity of parasites in each cell is less than 2. 

We now fix the growth rate of parasites g such that the probability that the Feller diffusion 
(X t : t > 0) reaches before it reaches 2 is strictly less than 1/2. Then the number of infected 
cells follows a supercritical branching process and grows exponentially with positive probabil- 
ity. Conditionally on this event, the proportion of infected cells does not tend to zero since 
the non-infected do not divide. Thus the organism does not recover. 

This example shows that we can have both non recovery and non proliferation of parasites in 
the cells: this corresponds to the 'moderate infection' we define in Section [522] 

Second, we focus on the linear division rate: 

r(x) — ax (x > 0), 

for some constant a > 0. Then a new cell appears at rate aX t = X^igvi r (^Q ) where we recall 
that X t , the total quantity of parasites at time t, follows a Feller diffusion with drift g > 0. 
As a consequence, as soon as the parasites do not become extinct, they grow exponentially 
and the number of cells grows with the same exponential rate. Thus the number of cells and 
the quantity of parasites grow similarly, which corresponds to non explosion of the average 
number of parasites per cell (see Section 15. 2. 3D . 

5.2.2 Proportion of cells moderately infected 

Using an approximation of the evolution of the proportion of cells based on Prop. 15 - H we can 
conjecture the following criterion for 'moderate infection'. We are interested in the case where 
the proportion of cells infected by more than A parasites vanishes (uniformly in time) as A 
tends to infinity. This means that the quantity of parasites in a cell chosen uniformly at time 
t does not explode as t — > oo. The criterion still depends on the maximum division rate 

r* = sup{r(a;) : x > 0} 

and on the sharing of parasites given by the random fraction 0. It gives an analogue of 
the recovery criterion given in the previous section for constant division rate. Again the bias 
phenomenon favors lineages with a large number of divisions, which explains the factor 2 in the 
criterion. Roughly speaking, if there exists a level of infection Xq beyond which the cells divide 
fast enough, then the parasites can not proliferate in a positive proportion of the cells and the 
infection is moderated. This level comes from the constant rate case: g < 2r(xo)E(log(l/0)). 
Simulations are provided as an illustration in the Section [5] 

Conjecture 5.1. We assume that r increases. 

(i) If g < 2r*E(log(l/0)) 7 then we have the following almost sure convergence: 

lim hm #HeV t :Xi>A} = ^ 

A->oo t— >+oo Nf 

(ii) If g > 2r*E(log(l/8)), then for every A > 0, 

{l im sup #{* 6 V t .Xl > A] >0} = { y t>QXt>0} as 

t->+oo N t 
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Let us here give some details on the approximation which leads to this conjecture. 
First, let us prove that almost everywhere on the set {limt^+oo X t = +00} we have limt_»+oo N t 
+00. Let A > 0, B = {lim t _> +00 X t = +oo;sup (>0 Nt < A} and x\ > such that r(x 1 ) > 0. 
Conditionally on the event B, there exists a random time r > such that Mt >r,X t > A.x\. 
Then for every t > r, there exists one cell with more than x\ parasites, and a new cell is 
created with rate at less r(x\). Hence, the number of cells is stochastically lower bounded by 
a Poisson process of rate r(x\). So P(B) = and this proves the result. 

When lim t ^ +oc Nt = +00, the bracket of the martingale part of (|5.2p converges to zero and 
the random fraction N t /(l + N t ) converges to 1. Thus we neglect them and consider the 
following deterministic approximation. For every / : ((t,x) 1— > ft(%)) £ Cl' 2 (R+ x *™ 



(f*t,ft)=fo(xo)+l / [ ^(s^) + gx—(s,x) + xa—^(s,x) ) fj, s (dx)ds 



J«, + 
H I q J r 2r(x) 

Jo 



ds 



dx 



f{s,0x)K(<to)-f(s,x) 



f(y) K (Vs,dy) - f(s,x) 



dx 2 
Jl a (dx) ds 

Jl s (dx) ds. 



(5.7) 



It describes the law of the following non-linear jump diffusion process (Ct ■ t > 0) which is our 
auxiliary process in this case. The process (Ct ■ t > 0) is a Feller diffusion with the following 
additional multiplicative jumps. At time t, conditionally on £ t = x, the process jumps to &x 
with rate 2r(x), and with rate E(r(( t )) to a state St defined by 



\S t G dx) 



r(x)P(( t e dx) 
K(r(Ct)) ' 



We can neglect the resampling term given by the jump St and we get a Feller diffusion with 
multiplicative jump with increasing rate 2r (.) as studied in Section r3.2l The extinction criterion 
of Proposition 13. 31 for this process gives the conjecture for the criterion for moderate infection. 



5.2.3 Average number of parasites per cell 

We now state some sufficient conditions under which the behavior of the average number of 
parasites per cell is known. 

Proposition 5.2. We assume that r increases. 

(i) If r is convex then: 

supE(AV7V t ) < 00. 
t>o 

(ii) If g > r* , where r* has been defined in h3.18\) . then: 

lim E(X t /N t ) = +00. 

One can note that in the first case r* = 00. 

Proof. First, notice that X t /N t = (^t,x). Applying (|5.2[) to f(x) = x and using that K is 
symmetric gives: 



-(fi s ,r)(fj, s ,x) 



ds + M, 



1,/ 



M, 



2,/ 



Using Jensen's inequality and the fact that r(.) and x 1— ¥ xr(x) are convex: 

E((^,r)< Ms ,z)) >E(r((n„,x))(n a ,x)) > r(E(( Ms , x)))E((fi s , x)). 



(5.8) 
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Thus, as N s /(1 + N s ) > 1/2: 

E((fi U x)) <(mo,x}+^ ( ff -ir(E((^.,x))))E(^ a ,x»ds. 

Since t !—>• E((/it,a;)) is upper bounded by the solution of the differential equation y' = 
(.9 — r (y)/2)y started at (/xqj^) which is bounded since there exists x > such that My > 
xo,9-r(y)/2 < 0. 

Under the assumptions of (ii), since (N s / (N s + l))(fx s , r)(^ s , x) < r* x), we obtain from 
(I5~a for f = 0: 

E((j*r t x))>(jMi,x)+ [ (g-r*)E((fi s ,x))ds, (5.9) 
Jo 

which gives the result since g — r* > 0. ■ 
5.3 Some results for decreasing division rate 

When the division rate decreases, very infected cells tend to become more infected whereas 
low infected infected tend to divide more and get rid of their parasites. We provide in Propo- 
sition [S~3l the following criterion: as soon as a healthy cell appears, the organism recovers, else, 
parasites proliferate in every cell. 

We are not interested here in constant division rate and we assume that r is decreasing 
and that: 

3xi > 0, r(xi) < r(0). (5.10) 
This means that non infected cells divide faster than other infected cells. 

Proposition 5.3. Conditionally on the event {3t > 0, 3 i G Vt : X\ — 0}, the organism 
recovers a.s. 

hm = 1 a.s. 

Conditionally on the complement event {Vt > 0, V i £ Vt : X\ > 0} ; for every A > 0, 

#{i : X\ <A}*H?Q a.s. 

Proof. First, if there exists t > and i G 3 such that X^ = 0, then iV t ° := #{ i G V t : X\ = 0} 
grows exponentially with rate r(0) after time to in the sense that 

< lim inf e^ r(0)t A^ t ° < lim sup e~ r(0) *7V t ° < +00 . 

Let us recall the notation iV t * = #{i G Vt : X\ > 0} that we used in previous proofs. We now 
prove that the number iV t * of infected cells grows geometrically with a lower rate that under 
(15.101) is at most: 

r= sup {r(xi)A + (f - A)r(0) +ln(p t -x)} < r(0), 

0<A<t<l 

where p t = P Kl (X f > 0) < 1 for t > is the survival probability of a Feller diffusion (without 
jump) at time t starting from x\. 

The set of infected cells at time t can be partitioned into subfamilies according to the first 
time A at which one of their ancestors' infection is less than x\. In lineages where the infection 
remains more than x\, the division rate is upper bounded by r(x\). 

Let us consider a family of cells for which A < t. Before this time the number of infected cells 
in this family grows at most with rate r(x\). After this time, the parasite population issued 
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from this cell becomes extinct at time t with probability at least 1 — pt-x since it starts with 
at most Xi parasites at time A. Moreover, after time A, conditionally on the survival of the 
parasites, the number of infected cells grows at most with rate r(0). 
Hence, for all xq > and t £ [0, 1], 

E* (iV7) < x E x (e r ^\ Pt . x .e r ^ t - x ^ < e\ (5.11) 

where E A denotes the expectation with respect to the r.v. A. This ensures that iV t * /N t goes 
to in L 1 and then in probability. Following the proof in the case of constant division rate 
gives the a.s. convergence. 

Second, since for every A > 0, 

inf F Xa (30 < £ < 1, Jf t = 0) > 0, 

0<x <A 

we can follow (|3.16p and conditionally on the event 

limsup#{i : X\ < A} > 1, 

t— ^oo 

there exists t > and i <G 3 such that X\ = a.s. This ends up the proof. ■ 

The following question arises: when does a non-infected cell appear ? We provide a suf- 
ficient conditions so that it happens a.s., which gives a condition for a.s. recovery depending 
on defined in (|3.20p . 

Corollary 5.4. If g < r*E(log(l/ min(0, 1 — 0)), then the organism recovers a.s. 

Proof. We consider the quantity Y of parasites in the following particular cell line: at each 
division, we choose the less infected daughter cell. If the mother cell is infected by x parasites, 
then the cell divides with rate r(x) and the quantity of parasites in each daughter cell is 
respectively Qx and (1 — 0)x, so that the less infected one has min(0, 1 — Q)x parasites. 
Thus the process Y follows a Feller diffusion with multiplicative jumps by min(0, 1 — 0) 
with rate r(.). By Proposition 13.41 this process becomes extinctcd a.s. under the assumption 
g < r*E(log(min(0, 1 — 0)). Thus there exist t > and i e 3 such that X\ = and the 
previous Proposition ensures that the organism recovers a.s. ■ 

6 Simulations 

In this section, we illustrate on simulations the criteria that are exposed in this work. For 
this, we fix a, K{d9) and r{x) and let the growth rate g of the parasites vary. For each values 
of this growth rate, we perform simulations and compute a relevant quantity (for instance the 
probability of extinction of the parasites) . 

Example 1: Let us first consider the case of a constant division rate r. We simulate the process 
(£t,i £ R+) that provides in this case an approximation for a lineage chosen uniformly. We 
consider a random fraction 0, uniformly distributed on [0, 1]. Then: 

E(log(l/0)) = / \og{l/x)dx = - [ \og(x)dx = -[x\og(x) - x]l = 1. 
Jo Jo 

In this case, Theorem 14.21 yields that if g < 2r the organism recovers a.s. Otherwise, parasites 
proliferate exponentially with positive probability. We see indeed on the simulations of Figure 
[T]that recovery happens with probability 1 when g < 2r and when g > 2r, the mean quantity 
of parasites grows exponentially. 
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Figure 1: Case of a constant division rate: r = 3, a = 2. eac/i division, the random fraction 
O is drawn uniformly in [0, 1]. We simulate for each value of g = b — d varying between 0.1 and 
10 (abscissa), N — 100 independent simulations of (£t,t £ [0,10 000]). Their values at T — 10 000 
provide Monte-Carlo approximations (a) of the extinction probability and (b) of the mean quantity of 
parasites on a log-scale. 
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Figure 2: Median (Thick plain line), mean (continuous line) and quantiles (2.5%, 25%, 75%, 97.5%, 
dotted and dashed) of the distribution of the proportion of cells infected by more than A — 10 parasites. 
The bound for the division rate is r* — 2 and K(d0) = 5\/^{dff). We have simulated the cells up to 
time T — 700,000. We let g — b — d vary between 1 and 3.7. For each value of g, we simulate 50 
branching Feller diffusions and compute for each simulation the proportion #{i € Vt : X\ > A}/N t . 
This provides an approximation of the distribution of this r.v. for the growth rate g. (a) r = r* = 2. 
(b) r(x) = r*(l - exp(-x/10)). (c) r(x) = r*l [5j+00 )(x). 



Example 2: Now, let us illustrate the Conjecture 15.11 that has been made for the case of a 
variable increasing division rate r(x). We consider the distribution of the proportion of cells 
infected by more than A parasites, #{« £ Vt : X\ > A}/N t , for constant and variable division 
rates. The real branching diffusions are simulated. Again, we have fixed a, K(dd) and r(x) 
and let g vary. 

We represent in Figure [5] (a) the case of a constant division rate. In the simulations of Figures 
[5](b) and (c), the rate is variable. In these cases, when the infection is low, so is the division 
rate. The infection hence lasts with a probability that is higher than in the case of Figure [5] 
(a). 

We can see different behavior depending on whether g is smaller or larger than 2r*E(log(l/0)), 
particularly when considering the evolution of the 97% quantilc. 

Although the use of the auxiliary process (ft : t > 0) fSection l5.2.2p for proving Conjecture 
15.11 has not been fully established yet, the role of the threshold 2r*E(log(l/0)) can be seen 
on simulations of f (Figure [3] (a)). Depending on the two regions, we see two regimes for the 
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Figure 3: Case of a non constant division rate: r(x) = 32(2,+aa)( x ): cr = 2. At each division, the 
random fraction is drawn uniformly in [0, 1]. We simulate for each value of g = b— d varying between 
0.1 and 10 (abscissa), N — 1000 interacting particles whose empirical distribution approximates the 
law of (Ct,« £ [0, 10 000]) (see the algorithm in JBGTlOf ). Their values at time T = 10 000 provide 
(a) extinction probability, (b) mean quantity of parasites on a log-scale. 



extinction probability. 



7 Appendix: Proof of Proposition 12.21 



Proof. We follow Fournier and Meleard [FM04j and break the proof into several steps. Recall 
that Z n is the microscopic process described in Section |2~21 and that we are under the assump- 
tions of Prop. 12.21 In the sequel, the constants C may change from line to line. 

Preliminaries: We recall moment estimates that will be useful in the sequel. Recall that Z n is 
started with one cell containing [nxo]/n parasites and that we assumed that r(x) < f(l + x p ). 
For all T € R+, using computation similar to |FM04| and stochastic calculus, we prove that 
for q eW: 

sup sup E((Z™,1 + x q )) < +00 and then sup E( sup (Z™, 1 + x p )) < +00, (7.1) 
lies* te[o,T] neN* te[o,T] 

by using Doob's inequality and the first part of (|7.1|l with q = 2p — 1. Let us also notice that: 

(Z?J)-(ZZJ) (7.2) 



JS.+ J [0,1] 

+ 



r(x) [/(&c) + /((l - - f(x)] K(d8)Zl l (dx)ds + M^ nJ + M?' nJ 



\j(x + i) - /(or)) (na 2 + 6)nx + (f(x - i) - /(a:)) (n<r 2 + d)» 



c/.s. 



where Af 1 '-^ and M 2 '-' are two square integrable martingales starting from and with brackets: 



(M 2 ^) t 



(f{x+^)-f(x)j {na 2 + b)nx 



f(x--)-f(x)\ (na 2 + d) 



Zl l (dx) ds 



[0,1] 



r(x) [f(6x) + /((l - 6)x) - f(x)} 2 K(dO)Zf(dx) ds, 



(M 1 ' n ' / ,M 2 '"' / ) t = 
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We shall prove that (Z n : n G W) is tight in D(R+, M F where .M F (R+) is embedded 
with the topology of weak convergence by using a criterion due to |EK86j . We will then con- 
sider the uniqueness of the limiting values of (Z n : n 6 N*) by identifying them as solutions 
of a certain martingale problem. 

Step 1: Tightness of (Z n : n G N*) Let us establish tightness on D([0,T],X F (R + )). For this 
we apply Theorem 9.1 p 142 of Ethier and Kurtz EK86 . Wc begin to prove that this sequence 
is tight in E)([0,T],.A/(f(R+)) with the vague topology on Mp(M.+ ). For this, we remark 
that the set of functions / G C^(R+,R) is a dense subset of Cb(R+,R) in the topology of 
uniform convergence on compact sets and prove that for such functions /, the sequences 
({Z n ,f) : n G N*) are tight in B(R + ,R). In a second time, we prove that the following 
compact containment condition holds: VT > 0, V?] > 0, 3K v ,t compact subset of .M F (R+), 

inf P (Z? G K-t, for t G [0, T]) > 1 - n. (7.3) 

Let / G C^(R + ,R). Using the Taylor- Young formula, there exists u™(x) and u^(x) in [0,1] 
such that: 

lim [x + - f(x)j (no 2 + b)nx + (f fx - - f(x)j (no 2 + d)nx 
hm (6 - d)xf(x) + ^(f( x+ ^m) + f' (x - ^ 



n— >-j-oc 



.C 



U l( X )\ , A t» ( „ U 2( X ) 



+ ^ bf" X +^l)+ d f" X 



2n \ \ n J \ n 

= (b-d)xf'(x)+o 2 xf"(x). (7.4) 

Under the assumptions of Prop. 12.21 the finite variation part V n ^ of (Z n , /) satisfies: 

K' f \ <nf\\oofT sup (Z?,l+a?) 
se[o,x] 

+ (\\f'\\oo(b-d)+o 2 \\f"\\ 00 + b -±^\\f"\\ 00 )T sup (Z?,x). (7.5) 
V Zn J se ro,T] 



By CCD, 



Using that 



sup E sup \V nJ \ <+oo. (7.6) 

nSN* \i£[0,T] / 



2 



f(x+±)-f(x)) (no 2 + b)nx+ (f (x - ij - /(a;) J (no 2 + d)nx 

= 2o 2 xf' 2 (x) + C{nAd ^ J)x , (7.7) 
n 

where C(n, &, d, cr, /) = O(l) when n — s> +oo, we obtain in the same manner that: 

supE sup \(M 1 ' n ' f ) t + (M 2 - nJ ) t \ ) < C sup E sup (Zf, 1 + x p ) ) < +oo, (7.8) 
neN* y*e[o,T] J neN* Vte[o,T] / 

by (|7.1[) . Let <5 > and let ((£„, T„) : n G N*) be a sequence of couples of stopping times such 
that S n < T n < T and T n < + <5. In the same way that we proved (|7.5[) . we can show that: 

EflV^-V^l) <C(b,d,o,r,f)5 sup E sup <Z t ",l + ] < O(b,d,o,f, f,T)6, (7.9) 
V " " / n£N* \te[o,T] / 



25 



by (|7.ip . The upper bound can be as small as we wish with a proper choice for 5. Similarly: 

E (\(M 1 ' n ' J ) Tn - (M 1 > nJ } Sn + (M 2 ' n ' f ) Tn - (M 2 ' n 'f) Sn \) 

<0(b,d,a,f,f,T)5 sup E [ sup (Z?,\ + x p )\ 



nGF 



. tG[0,T] 



<C(b,d,a,f,f,T)S. 



(7.10) 



Then Aldous-Rebolledo and Roelly's criteria |JM8 6, Roe8Sj ensure that the sequence (Z n ) ne ^ 
is tight in D(R + , A4f(R+)) where .M F (R + ) is embedded with the vague convergence topology. 

Let us now prove the compact containment condition (|7.3|) to obtain the tightness in 
B(R+, Aip (&.+)) with the weak convergence topology on .M F (R+)- Recall that the sets 
M<N o ([0,ao]) °f measures with mass bounded by Nq and support included in [0,ao] are 
compact (see |Kal83| . Sect. 15). Notice that: 

{Z? i M< No ([0,a a }), t G [0,T]} c {3t G [0,T], N? > N Q } U {3t G [0,T], X t n > a } 

where iV" is the number of cells and where X" is the total quantity of parasites at time t. 
Hence: 



(Z? i M< No ([0,a o }), t G [0,T]) <P(3i G [0,T], iV™ > A ) +P(3t e [0,T], X? > a ) 



<_L E ( sup N r t 



'Nn 



te[o,T] 



1 

a 



E( sup X\ 

te[o,T] 



(7.11) 



Thanks to (|7.ip for fixed T, we obtain (|7.3p by choosing iVo and ao sufficiently large. This 
concludes the proof of the tightness of {Z n : n e W) in D([0,T],X F (M + )). 

Step 2: Identification of the limit Let us consider an adherence value Z of the sequence (Z n : 
n G N*), and let us denote again by (Z n : n G N*) the subsequence that converges towards Z 
in law in D([0,T],Mjr(l + )). Let / G (R+,R). For k e N*, let < i x < • • • t k < s < t < T 
and <pi,--- ,(p k e C b (A4 F (R+),R). For z G D([0,T],X F (R+)), we define: 



= Plfoi)' ' • <Pk(z* h ) (Z t ,f) - (z s ,/) - 



(r(s) (/(&:) + -/(*)) 



'R+ «/[0,X] 

+ xf'{x)(b -d) + <7 2 xf"(x))K(d6)z u (dx) du\ . (7.12) 

Then |E (*(Z))| < A + B + C where: 
A = \E(^(Z))-E(^(Z n ))\ 

E (*(£")) - E (MK) ■ ■ ■ MZ? k ) [Mt nJ ~ Ml^f + M^ nJ - M 2 ^] ) 
E (y x ■ ■ ■ Vk {Zl\) [M^ nJ - M]- n ^ + M 2 ' nJ - M 2 ' n ' f ^ ) 



B = 
C = 



The map z G B([0, T], A4f(R + )) is continuous as soon as (ti, •••£&, s,i) does not 

intersect a denumerable set of points of [0, T] where Z is not continuous (e.g. Billingsley 
Bil68] Theorem 15.1 p. 124). The convergence in distribution of Z n to Z, together with (|7.1|) . 
implies that A converges to when k — > +oo. Since M 1 '™'-' and M 2 > n >* are martingales, C = 0. 
From (|7~j]) . 



S < 



E 



-x 



x 
2ri 



+ 

bf"[ 



f" (x + 



u^{x) 



df" x 



+ f" [x 



u?(x) 



2/"(x) 



Z n s (dx)ds 



2G 



Since / is of class C 3 , the integrand is upper bounded by Cx/n. Using (|7.1[) , B < C'/n. This 
proves that K(^!(Z)) — and hence W* is a martingale. The computation of its bracket is 
standard and obtained by taking the limit in (Z™,f) 2 — (Zq,/) 2 on the one hand, and by 
using Ito's formula with (|2.5[) on the other hand. 



Step 3: Conclusion In the Step 2, we have identified the adherence values of the sequence of 
processes (Z n : n G N*) as the solutions Z of the martingale problem associated with the 
following generator A. Following Dawson |Daw93| . Section 6.1, we choose for the domain 
V(A) the set of functions of the form F f (Z) = F((Z,f)) with F G C 6 2 (M,IR), / G C fc 2 (M + ,]R) 
and where Z G Mf Let: 

AF f (Z) =F'((Z, /)) f [(b - d)xf(x) + <J 2 xf"(x)] Z(dx) (7.13) 
+F"((ZJ)) [ a 2 xf' 2 (x)Z(dx) 

+ f f r(x) \f((Z, J) + f(9x) + /((l - 6)x) - f(x)) - F f (Z)] K(d6)Z(dx). 

JR+ J [0,1] 

The generator A is linear. Since it is the infinitesimal generator of the process (|2.3[) - (|2.5p . A 
is closed by Corollary 1.6 p. 10 of |EK86| . This implies that the resolvent set of A contains 
(0,+co). By Proposition 3.5 p. 178 of |EK86j . since for any initial condition Zq G .Mf(R+), 
(|2.3[) - (|2.5[) defines a solution of the martingale problem, the generator A is dissipative. These 
last facts imply by Corollary 4.4 p.187 of [EK86] that two processes of D(K + , X F (M+)) 
satisfying the martingale problem associated with A have the same distribution. ■ 
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